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ABSTRACT 

We present a pair of 3-d magnetohydrodynamical simulations of intermittent jets from a central 
active galactic nucleus (AGN) in a galaxy cluster extracted from a high resolution cosmological sim- 
ulation. The selected cluster was chosen as an apparently relatively relaxed system, not having un- 
dergone a major merger in almost 7 Gyr. Despite this characterization and history, the intra-cluster 
medium (ICM) contains quite active "weather" . We explore the effects of this ICM weather on the 
morphological evolution of the AGN jets and lobes. The orientation of the jets is different in the two 
simulations so that they probe different aspects of the ICM structure and dynamics. We find that 
even for this cluster that can be characterized as relaxed by an observational standard, the large-scale, 
bulk ICM motions can significantly distort the jets and lobes. Synthetic X-ray observations of the 
simulations show that the jets produce complex cavity systems, while synthetic radio observations re- 
veal bending of the jets and lobes similar to wide-angle tail (WAT) radio sources. The jets are cycled 
on and off with a 26 Myr period using a 50% duty cycle. This leads to morphological features similar 
to those in "double-double" radio galaxies. While the jet and ICM magnetic fields are generally too 
weak in the simulations to play a major role in the dynamics. Maxwell stresses can still become locally 
significant. 

Subject headings: galaxies: jets - galaxies: clusters: general - methods: numerical - X-rays: galaxies: 
clusters - magnetohydrodynamics (MHD) 



1. INTRODUCTION 

Observations of powerful outflows from active galac- 
tic nuclei (AGN) show that they have a considerable 
effect on the energy budget and morphology of their 
host galaxy cluster. The jets from the super- massive 
black hole that powers an AGN inflate low density bub- 
bles in the host galaxy cluster's intra-cluster medium 
(ICM) that are observed as X-ray cavities. The mini- 
mum energy required to inflate these cavities can be as 
high as 10^^ erg {e.g., McNamara et al. 2005 : Wise et al\ 
[2007) and typically imply a power on the order of 
10^^ - 1 0^^ erg s~^ a veraged ove r the age of the cav- 
ity (e.(/.,[Birzan et al. 2004; McN amara fc NulsenI [20071 : 
iBirzan et al\ 120081 : [olehl et aL 120081 ). Most X-ray cav- 
ities are paired as a result of the twin AGN jets that 
formed them, and most are filled with radio synchrotron 
emission from relativistic, cosmic ray electrons (CRs); 
that is, "radio lobes." Absent thermal conduction, mod- 
els of the ICM without any heating source predict that 
galaxy clusters should establish "cooling- flows" , where 
radiatively cooled ICM plasma falls towards the central 
galaxy drivin g up star formation to very high values 
(jFabianl 119941) . Observations do not show evidence of 
these classical cooling- flows, but in contrast they show 
that mass accretion onto the central galaxy is at least 
an order of m agnitude below cooling-flow predictions 
(jPeterson fc Fabian 2006 ). One popular explanation for 
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regulating cooling in clusters is that energy output from 
the AGN in the central galaxy heats the ICM sufficiently 
to limit accretion. Indeed, the enthalpy content esti- 
mated in X-ray cavities appears to be sufficient to offset 
cooling of the ICM {e.g., Birzan et al. 2004, 2008). While 
the details of the physical processes involved in the trans- 
fer of energy from AGN outflows to the ICM are not 
fully understood, it is generally agreed that these out- 
flows contribute to the dynamics and thermodynamics 
of the ICM. 

At the same time there is also good evidence that the 
ICMs significantly affect AGN outflows. For example, 
several radio lobes that have been ob served, such as Hy - 
dra A (jWise et a/.ll2007r ) and 3C 31 (|Laing e^ a/.l 120081 ). 
show complex morphologies with various kinks and turns. 
The cited authors suggested that, for these objects, the 
jets and lobes may be reacting to bulk motions of the 
IC1V0. Further, there are groups of AGN-produced ra- 
dio sources specifically classified by their bent morpholo- 
gies; namely, narrow-angle tails (NATs) and wide-angle 
tails (WATs), dependin g on th e apparent degree of bend- 
ing. lOwen fc Rudnicfl (|1976f ) proposed that these ob- 
jects obtain their bent shape from the relative motion 
through the ICM of the galaxy hosting the AGN engine. 
In the case of WATs, the low peculiar velocity of the 
associated galaxies in combination with reasonable as- 
sumptions for the density and velocity of the plasma 
within the AGN outflows suggest that the motion of 
the ICM within the cluster is more responsible for their 
characteristic C-shape than the motio n of the individual 
host galaxy within the cluster {e.g., iBurns et a/.l 11994] : 

^ In this work we use "bulk flow" or "wind" to mean motions of 
the ICM that are not characterized by isotropic turbulence. The 
complete motion of the ICM is described by the combination of 
bulk flow and turbulence on scales smaller than the cluster core. 



iDoudass et"a/]l2QQ8[ ). iBurns et al\ (|1994[ ) argued that 
cluster mergers can produce sufficient ram pressure in 
the ICM to cause this shape, and simulations have shown 
that motions produced by merger events with sufficient 
ram pressure to bend WA Ts are persistent for timescales 
exceeding 10^ years (Loken e^ a/.l 119951: iRoettiger et al\ 
119961 : iRicker fc Sarazinll2QQlf ). Similarly, recent simula- 
tions show large-scale ICM "sloshing" several gigayears 
after a gravitational encou nter with a passing subcluster 
(e.g., iZuhone et a/J|2010D. 

Direct detection of such bulk ICM flows became fea- 
sible with the Advanced Satellite for Cosmology and 
Astrophysics (ASCA), and the Chandra X-ray Ob- 
servatory has further im proved such measurements. 
iDupke fc BregmanI (|2QQ6l ). for example, used Chandra 
observations to explore the velocity structures in the 
Centaurus Cluster (Abell 3526), and they found that 
significant departures from hydrostatic equilibrium are 
possible even in a cool cluster that shows little evidence 
of a strong merger. ICM sloshing motions have also been 
observed in the central regions of seve ral cool-core clus- 
ters (see 'Markevitch fc Vikhlininll2QQ7l for a review). 

There is limited simulation literature addressing inter- 
actions between AGN outflows and dynamical ICM flows, 
or " weather . " Fo r exam ple , .Balsara fc Norman"! ([l992, ) 
and iPorter et all (|2QQ9[ ) carried out 3d simulations of 
AGN outflows into supers onic cro ss-winds that produced 
NAT-like structures. Loke n et all (1995) used 3-d hydro- 
dynamic simulations of a low density jet obliquely cross- 
ing a plane shock to demonstrate that merger produced 
shocks can decollimate a jet and cause it t o bend. In a 
study of AG N energy depos i tion to the ICM iHeinz at~^ 
|2006) and iMorsonv etal\ (|2QlQr ) used output from a 
cosmological simulation as initial conditions for hydro- 
dynamical simulations that included a pair of AGN jets. 
To differentiate behaviors fully from hydrostatic clusters 
the latter study specifically selected a cluster that was far 
from being relaxed and that included a net rotation and 
a dynamically induced cold front. Not surprisingly those 
authors found significant effects on jet propagation and 
lobe morphology due to the rather severe ICM weather. 

Especially in the context of clusters exhibiting cavi- 
ties produced by central dominant galaxies, which more 
typically associate with apparently more relaxed clus- 
ters, however, it is important to understand if such ef- 
fects can ta ke place in those types of clusters. Except 
for the Port er et all (2009) idealized NAT study, there 
are no published simulations of AGN interactions with 
dynamical ICMs that included magnetic fields. How- 
ever, since magnetic fields both introduce potential dy- 
namical influences and are also essential to modeling the 
generation and transport of diagnostic nonthermal emis- 
sions, it is imp ortant to extend s uch si mulations into the 
MHD regime. lO^Neih fc JonesI ([Moh (hereafter OJIO) 
recently conducted a series of 3-d MHD simulations of 
steady and intermittent jets in an analytically defined, 
magnetized ICM close to hydrostatic equilibrium. The 
magnetic field was disordered on scales similar to the 
size of the ICM core and on average produced a pres- 
sure roughly 1% of the gas pressure, so probably char- 
acteristic of cluster strength fields. They found that for 
localized regions, those magnetic fields can have a signif- 
icant effect on the morphology of the jet lobes. Whether 
this remains true for a dynamically self-consistent cluster 



evolved through a cosmological simulation has not been 
examined. Another important result from the OJIO sim- 
ulations was the lack of any significant deflections of the 
jets or lobes. The steady and intermittent jets in their 
study both maintained relatively straight-line paths, sug- 
gesting that a dynamical ICM is required to produce bent 
jet and/or bent lobe morphologies. 

To address some of these open issues we present here a 
pair of 3-d simulations of MHD jets in a dynamic, mag- 
netized ICM extracted from a high resolution cosmolog- 
ical simulation. The selected cluster was chosen from 
among the more relaxed systems in the simulation. To 
expand the insights into these interactions and the influ- 
ence of ICM weather on radio lobe morphology the two 
simulations directed the AGN jet pairs in two mutually 
orthogonal directions that probed distinct regions of the 
ICM. The paper is organized as follows: ^describes the 
numerical setup for the calculations, ^S] describes the ini- 
tial conditions in the cluster, g] qualitatively describes 
the results of the simulations, ^ discusses the effects of 
weather on the jets and lobes, and ^ summarizes the 
major conclusions. We assume below that Hq = 70 km 
s-^ Mpc-\ Qm = 0.3, and Qa = 0.7. 

2. NUMERICAL DETAILS 

The simulations of AGN jets for this study were per- 
formed on a 1008x1008x1008 three dimensional Carte- 
sian grid with a uniform resolution of Ax = Ay = 
Az = 1 kpc per zone using WOMBAT (Mendvgr al et all 
l2012af ). which is a versatile, efficient and highly scal- 
able MHD code. Here, WOMBAT evolved the equa- 
tions for non-relativistic ideal MHD with adiabatic in- 
dex, 7 = 5/3, using the second orde r, Total Variation Di - 
minishing (MHDTVD) method of Rvu fc JonesI (|1995aD . 
The solenoidal constraint of the magnetic field was main- 
tained to machine accur acy with the ^^con strained trans- 
port" (CT) method of iRvu et al\ (|1998[ ). which solves 
the induction equation with a directionally unsplit up- 
date using EMFs derived from the Riemann solver. The 
simulations included the total energy conserving method 
for gra vity available in WOMBA T and described in de- 
tail in iMendvgral et al\ (|2012br ). This method adds 
gravitational source terms to MHDTVD in a manner 
that maintains 2"^ order code accuracy. The evolu- 
tion of a passive cosmic ray electron (CR) momentum 
distribution, /(x,p, t), was also calculated in the sim- 
ulations with a restricted form of t he ^^Coa rse Grained 
Momentum Volume" (CGMV) (Jon es fc Kan g 2005) im- 
plementation in WOMBAT. CGMV solves the diffusion- 
convection equation for a population of relativistic par- 
ticles. For these simulations spatial CR diffusion out- 
side shocks was neglected, while adiabatic changes in 
momentum, test-particle diffusive acceleration at shocks 
and radiative losses by synchrotron emission and inverse 
Compton scattering of Cosmic Microwave Background 
(CMB) p hotons were included (see further details below 
and, e.g. JTredllis et al\ l\2001\ )). 

The boundary conditions for these simulations were 
based on a technique used by OJIO. To limit any in- 
fluence from the boundaries, we used a very large do- 
main (over a Mpc along a side) and modified continuous 
boundaries that Hmit the fluxes through those surfaces. 
In the boundary zones a hydrostatic atmosphere with a 
constant sound speed was derived by extrapolating con- 



ditions from the last physical zone. Since the mass distri- 
bution was known on a much larger domain (see below), 
the gravitational potential was known out to radii far 
beyond the computational domain. Consequently, the 
exact value for the potential described below was used 
through the boundaries. 

2.1. Setting up the Cluster 

The initial conditions for the simulations described 
here were extracted from a high resolution simulation 
of a galaxy cluster performed with an MHD implemen- 
tation of the smo othed particle hydrodyna mic (SPH) 
code GADGET-3 (iDolag fc Stasvszvri [2QQ9h The se- 
lected cluster is q676 from iDolag et al\ (|2QQ9l ), resim- 
ulated at extreme high resolution with MHD enabled 
(for details see Stasyszyn & Dolag, in prep). In short, 
the cluster used for this study has been extracted from 
a re-simulation of a Lagrangian region selected from 
a cosmological , low er resolution DM-only simulation 
(|Yoshida et al\ l2QQlf ). This parent simulation has a 
box-size of 684Mpc and assumed a baryon fraction of 
/bar = 0.13 and ag = 0.9 for the normalization of the 
power spectrum. Using the ^^Zoomed I nitial Conditions" 
(ZIC) technique (JTormen et al\ [TOQ?! ), this region was 
re-simulated with higher mass and force resolution by 
populating the Lagrangian volume with a larger num- 
ber of particles, while appropriately adding additional 
high-frequency modes drawn from the same power spec- 
trum. The initial unperturbed particle distribution (be- 
fore imprinting the Zeldovich displacements) w as real- 
ized through a relaxed glass-like configuration (|Whitel 
|1996[ ). Gas was then added to the high-resolution re- 
gions by splitting each parent particle into a gas and a 
DM particle. The gas and the DM particles were dis- 
placed by half the original mean inter-particle distance, 
such that the centre-of-mass and the momentum of the 
original particle are conserved. The final mass resolution 
of the dark matter and gas particles in our simulations 
is mdm = 12.1 X 10^ M0 and mgas = 2.2 x 10^ M© re- 
spectively. Thus, the clusters within its virial radius is 
resolved with 11 x 10^ dark matter and a similar num- 
ber of gas particles. In the simulations, the gravitational 
softening length was chosen to be e = 1.4 kpc which cor- 
respond roughly to the mean particle separation in the 
center of the cluster at z=0. For the initial magnetic field 
we choose a space filling, homogeneous, prin iordial mag- 
netic field oi B = 10~^^G co-moving as in iDolag et al\ 
(|2009h . 

For this work, we were interested in assessing the ef- 
fects of cluster "weather" specifically for a cluster that 
observationally appears relaxed yet harbors significant 
bulk ICM flows. We selected the state of g676 at a red- 
shift of z ~ 0, which was the least dynamic period for this 
cluster. As seen in Figure [H the last major merger (mass 
ratio of 0.27) occurred at z ^ 0.8 (a little less than 7 Gyr 
before our jet simulations start), and the dark matter 
mass continued to grow slowly after that. The general 
morphology and environ ment of th e cluster at z = can 
be seen in Figure 1 of Dol ag et al\ (I2QQ 9). Volume aver- 
aged values for baryonic density, temperature, velocity, 
and the magnetic field were all computed from the SPH 
particles and mapped to a uniform grid matching the 
jet simulation grid. The simulation assumed a hydrogen 
fraction of X = 0.76 and a helium fraction of F = 0.24, 



resulting in a mean molecular weight of /i = 0.59. An 
ideal equation of state was used to convert between pres- 
sure and temperature. Radiative cooling of the thermal 
plasma was not included in the GADGET calculations 
that formed cluster g676^ so neither was it used in the jet 
simulations. A summary of relevant cluster ICM proper- 
ties is given in ^ 

For the jet simulations, we used a static gravitational 
potential defined from the total gravitating mass (dark 
matter and baryonic matter) azimuthally averaged from 
the SPH values at the time used for the initial condi- 
tions. Figure [2] shows the gravitating mass as a function 
of radius. The total enclosed mass out to the virial ra- 
dius of 1.4 Mpc is 1.53 X lO^'^M©. The decision to use 
a static potential as opposed to evolving a population of 
DM particles and maintaining an updated full potential 
was based on simplicity, the relaxed nature of the clus- 
ter and the relatively brief duration of the jet activity 
(200 Myr) compared to the dynamical timescale for the 
cluster mass distribution. Obviously over longer times 
or in a cluster involved in significant active restructur- 
ing it would be important to follow the evolution of the 
gravitational potential. 

2.2. Setting up the Jets 

For each simulation, a pair of oppositely directed jets 
was launched from a cylinder centered in the grid (X = 
y = Z = 0). The launching cylinder had a radius of 
Vj — 3 kpc and length on each side of the jet origin of 
Ij = 8 kpc. An additional 3 kpc thick collar and cap 
surrounded the cylinder to transition from jet conditions 
to ICM conditions and to contain the return electric cur- 
rent associated with the jet magnetic field ( g2.2.2p . Spe- 
cial consideration was given to how the cylinder inter- 
acted with the ICM. In these simulations, ambient fiows 
could in principle enter the cylinder, even from the initial 
conditions. This can be particularly problematic if ICM 
magnetic fields enter the cylinder, where they can blend 
with the jet field and produce undesirable field geome- 
tries. 

Two steps were taken to alleviate this. First, at the 
start of the simulation the ICM material in the vicinity 
of the jet cylinder location was gently pushed aside be- 
fore the jet cylinder was initialized. In addition, refiect- 
ing boundaries were defined as defaults for the cylinder. 
The initial ICM evacuation was accomplished through a 
purely solenoidal velocity field imposed at the start of 
the simulation on a sphere with radius Vj centered on 
the grid with v^ = Ve{r/rj) , and Ve set to the local 
sound speed. This sphere was over-pressured by 90% as 
well. These conditions were allowed to evolve for 3.28 
Myr, which was the sound crossing time over Tj. This 
removed ^ 25% of the ICM plasma and magnetic energy 
from the jet cylinder region. The jet cylinder was then 
established at t = 3.28 Myr. To divorce the ICM and jet 
magnetic fields, any remaining ICM field penetrating the 
jet cylinder or collar was set to zero at this time. This 
last step effectively introduced some magnetic monopoles 
around the cylinder, but they had negligible effect on the 
jet fiows and were maintained to low levels by the CT 
scheme in WOMBAT. The collar region around the ac- 
tive jet cylinder acted as an impenetrable electrically con- 
ducting surface because of the imposed refiecting bound- 
aries for both the normal momentum and magnetic field. 
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Figure 1. Evolution of velocity dispersion (dark matter) and total mass as a function of redshift for the selected galaxy cluster g676. The 
vertical lines mark identified merger events with a mass ratio indicated by the labels. 
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Figure 2. Gravitating mass profile used at all times to define the gravitational potential. The mass includes both dark matter and 
baryonic matter from the GADGET simulation and was derived from an azimuthal average at each radius. 



This prevented loss of any conserved quantity into the 
launching region. The caps of the cylinder transitioned 
between reflecting boundaries and outflowing jet proper- 
ties as the jet cylinder transitioned between active and 
inactive states. 

An important feature of the jet cylinder for this inves- 
tigation was that it could be oriented in any direction 
on the grid. This allowed jets to be launched at an ar- 
bitrary angle with respect to the large scale flows and 
pressure structures in the ICM. To achieve this, we de- 
flned all vector quantities in the jet cylinder according to 
cylindrical coordinates. The jet cylinder is azimuthally 
symmetric, so only two coordinates are needed, z and r, 
where z is along the cylinder symmetry axis and r is radi- 
ally outward. By default the jet z axis was aligned with 
the grid Z axis, and Euler angles were used to reorient 
the cylinder. After the rotation, all cylinder coordinates 
were converted into grid coordinates so that quantities 
could be mapped to the Cartesian grid. The speciflc ori- 
entations of the two jets in this study are discussed in 
section ID 

The active jet density and pressure were set to pj = 
4.165 X 10-2^ g cm-3 and P^ = 2.499 x 10"^° dyne cm-^ 
(approximately 1/200 the initial local ICM density) and 
in approximate pressure equilibrium with the initial sur- 
rounding ICM. For the assumed adiabatic index of 7 
5/3, this gave a sound speed inside the jet of c^ 



small contributions from magnetic energy, is given by 



- 10^ 



km s~^ or 0.03c. A passively advected "color" (or mass 
fraction) tracer, Cj, was also introduced in the jet cylin- 
der that was used to identify AGN and ICM plasma. 
The color was set to 1 in the jet cylinder, while all ICM 
material was initially assigned Cj =0. 



2.2.1. Jet Velocities 

When it was active the jet velocity in the jet cylinder 
was given a spatial ramp from its origin in z to allow a 
smooth transition between the oppositely directed flows 
of the two jets. In particular, the velocity had the form. 



j{z) = VjT{z)^ where 



T{z) = TANK 1.^ 



(1) 



The value 1.8 was determined empirically by minimizing 
the automatic apphcation of diffusive, positive-deflnite 
HLL fluxes by WOMBAT in the jet cylinder, which are 
used to eliminate unphysical pressures during peak flow 
times. The velocity in the jet cylinder was deflned with 
this ramp to be Vzj{z^t) = Vj{t)F{z)^ where Vj{t) = 
1.2 X 10^ km s~^, when the jet was fully on (see ^2.2.3p . 
This yields a pair of supersonic emergent jets with the 
flow at Mach number, Mj = 1.2, with respect to the 
jet sound speed. The peak speed of Vj corresponds to a 
Lorentz factor of 1.0008 and is sufficiently low that a non- 
relativistic MHD solver is applicable. The jet velocity 
was uniform in r across the jet core (r < 3 kpc), with a 
transition to the ambient medium velocity ex r~^/^. 
The core luminosity of each of the two jets, ignoring 



Lj =7rrjVj -pjv 



1 



7 



2^'^ ' 7-1 



= Trr.-pV^ 1 



Mf{j-i)r 



(2) 



where 7 = 5/3. When the jets are fully active, the pre- 
ceding parameters give a total luminosity (including both 
jets) of Lj = 6 X 10"^^ erg s~^. The transitional collar 
region adds roughly another 50% to the actual jet lumi- 
nosity. The ratio of kinetic energy ffux to enthalpy ffux 
for these jets was pjv'j/bPj = 0.48. The actual total en- 
ergy injected onto the grid by the jet cylinder in these 
runs was 2.2 x 10^^ erg. 

For later discussion it is also useful to know the thrust 
of the jets, Fj. The thrust is easily computed from the 
luminosity (with 7 = 5/3) as 



Fj 
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SO that the individual core jet thrusts at peak luminosity 
are Fj ^ 1.6 x 10^^ dyne. 

2.2.2. Jet Magnetic Field 

The magnetic field in the jet launching cylinder was 
derived with a method that maintained V • B to machine 
accuracy. The jet field was purely toroidal inside the jet 
cylinder with 5^ = Bjir/vj) when r < r^, where Bj = 
7.92 pG (plasma p = SttPj/B'^ = 100). Outside of that 

radius, the field dropped off as r~^ and was zero outside 
of the jet cylinder collar. This kept the jet field separated 
from the ICM field and made the net electric current in 
each jet zero. The arbitrary orientation of the jet cylinder 
was an additional complication for defining the jet field. 
We implemented this by defining an electric field that 
determined the rate of change of the magnetic field in 
the cylinder according to Faraday's Law, —d'B/dt = V x 
E. Solving this equation allowed magnetic field to be 
replaced as it was advected out of a zone at an arbitrary 
angle. Given the desired form of the magnetic field and 
the jet velocity as a function z, this electric field was 

'Bi-.-W2fe{(77)'-[f-K^)']} 

for r < rj and z\ < lj 

for rj<r< r^ and \z\ < lj 



E,{r,z) 







otherwise. 



(4) 



where r^ = Tj + 3Az. The magnetic field added to the 
zones in the jet cylinder and collar was simply ^B = 
At(V X E), where At is the current time step. The value 
of V X E only needed to be computed once at the begin- 
ning of the simulation. 

2.2.3. Jet Source Time Evolution 

AGN activity is likely to be unsteady on the > 10^ yr 
timescales required to form large scale cavities, while the 
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intermittency of jet activity has significant impact on 
the interactions with the ICM (see, e.g., OJIO). In these 
simulations the activity of the jets was toggled on and 
off over approximately 26 Myr as a simple way to mimic 
such cycHc activity. In particular we applied to the char- 
acteristic jet velocty, Vj, and magnetic field, Bj, a ramp 
function, R{t)^ defined as 



R{t) = r]{t)MAX 

{l-r]{t))MAX 
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(5) 



(6) 



where ton and toff are the start and end time for the 
current jet cycle and tbiend = 3.28 Myr. A duty cycle 
of ~ 50% was used for both simulations with a period 
of 26.2 Myr. Thus, the jet velocity and added magnetic 
field were given a time dependence as Vj (z, t) = Vj {z)R{t) 
and 5B{r,z,t) = At[V x E{r,z)]R{t). 

2.3. Cosmic Ray Electrons 

In order to enable realistic synthetic observations of 
nonthermal emissions we incorporated in the jet outflows 
a population of test particle CR electrons with Lorentz 
factors from 10 to 1.6x10^ and evolved their distribution, 
/(r,p,t), using the CGMV routine in WOMBAT. This 
range in Lorentz factors is sufficiently large for the typi- 
cal magnetic field strengths in these simulations to model 
synchrotron emission from a few MHz to tens of GHz. No 
CR electrons were included in the initial ICM, although 
any electrons introduced by the jets and then mixed into 
the ICM were followed in the CGMV scheme. For these 
simulations, CRs were subject to adiabatic changes in 
momentum as well as energy losses due to synchrotron 
emission and inverse Compton scattering of CMB pho- 
tons. CRs were also subject to 1®* order Fermi accelera- 
tion at shocks (diffusive s hock accelerat ion (DSA)). 

As described by Tr egillis et al.\ (|2001 '). the character- 
istic scattering length for electrons in this energy range 
and in field strengths of a few /iG, typical for these simu- 
lations, is on the scale of the solar system, and the accel- 
eration timescale is on the order of years. These scales 
are both unresolved in these simulations, which resolves 
lengths of a kiloparsec and time in kiloyears. This mis- 
match in scales can be exploited by assuming that the 
acceleration process for CRs at a shock is completed in a 
time that is far less then a typical dynamical time step in 
the simulation. For this application, the CGMV imple- 
mentation does not include diffusion terms and simply 
redistributes existing post-shock CRs into a power law 
according to test particle DSA theory. The slope Qs of 
the resulting power law is related to the measured com- 
pression ratio of the shock, r, by qg = ;^. To implement 
DSA WOMBAT includes a directionally unsplit shock 
detector that locates and characterizes shock strengths. 
Eight boundary zones along domain boundaries allowed 
accurate shock detection and characterization up to the 
edge of the local domain boundary and the physical grid 
boundary. 

The CR population was introduced in the flow of AGN 
plasma out the jet cylinder. Since the CR electrons were 



passive, their number was somewhat arbitrary, except in 
the calculation of nonthermal emissions. For the syn- 
thetic observations reported here, we assume the emer- 
gent jets carried a population of CRs with a density of 
5 X 10~^ cm~^ (about 0.2% of the jet bulk proton num- 
ber density) and had a single power law slope of qj = 
4.5. For consistency this population should not be ex- 
pected to contribute significantly to the subsequent flow 
dynamics. In fact, with this density and slope, the nom- 
inal CR electron pressure at the jet source was 4 x 10~^^ 
dyne cm~^, or approximately 2% of the "thermal" gas 
pressure in the AGN plasma. Especially, given the softer 
equation of state for these CRs, the downstream CR pres- 
sure influences would, indeed, have been minor. Since 
the nonthermal emissions of interest to us here are as- 
sociated with the jets, we disabled injection of fresh CR 
electrons from the thermal pool for these simulations in 
order to avoid contamination by CRs that would have 
been injected in ICM shocks. 

2.4. Synthetic Observations 

Synthetic observations of the simulations were calcu- 
lated following t e chniq ues similar to those described in 
iMendvgral et all ()2011l ). We focus on two models for op- 
tically thin emission; thermal bremsstrahlung from the 
hot plasma and synchrotron radiation from the CR elec- 
trons. X-ray emission in an energy band typified by 
Chandra and low frequency radio emission similar to 
those obtainable by LOFAR are of particular interest. 
We note that, although emission due to inverse Comp- 
ton scattering of CMB photons is simple to include, it is 
negligible for the energy/frequency range of interest with 
the CR electron numbers we assume here. 

The thermal bremsstrahlung emissivity was calculated 
in every computational zone as 



5.4 X 10 ^^gff{iyiocahTe)x 



(7) 
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where viocai = ^obs (1 + ^) for a redshift of z, Tg is in keV 
and all other quantities are in cgs units. The free- free 
Gaunt factor, gff^ was computed by interpolation from 
the values calculat ed for plasma w ith typical ICM prop- 
erties in Table 1 of lNozawa et al\ ([T998). For simplicity, 
we assumed a fully ionized, Zi = 1 pure hydrogen gas 
with an ideal equation of state. The pure hydrogen as- 
sumption, however, does not match the abundances that 
were used to d efine the composition of the baryons de- 
scribed in §2.1[ Although formally the Gaunt factor does 
depend on the ion charge Z^, the modifications to gjf for 
typical temperatures and observations energies of these 
simulations for a Zj = 2 g as were only a few percent from 
Table 1 of Nozawa et al\ (p^998) . The remaining depen- 
dence is the Zf term in Equation [8l Since the abundance 
is uniform on the grid, this only sets a normalization 
for bremsstrahlung emission. Finally, only for the ther- 
mal emission computation, the mean molecular weight /i 
= 0.5 was used to define the temperature in a zone as 
Te = Ti = T(keV) = /iPm///(1.602 x lO'^p), where P, 
rriH and p are all in cgs units. To simulate the band- 
width of real instruments. Equation [8] was numerically 
integrated over a range of frequencies. Unlike the tech- 
nique used in ^Mendygral et a l., ( 2011. ) , where emission 



from zones with AGN plasma was clipped for numerical 
reasons, we do not eliminate emission from those zones 
in this study. Significant and complex mixing between 
ICM and AGN plasma in these simulations prevented a 
clear cutoff based on the passive color tracer, Cj. 

Synchrotron emissivity was c alculated us i ngs; the local 
CR spectrum and was given bv I Jones et al\ ()1974l ) as 






,cx+l 



local 



(8) 



where /p^ is the CR phase space density at momentum 
Pc, (Ic is the CR spectral slope at pc and jc is a constant 
of order unity. The spectral index is related to the CR 
spectral slope as a = {qc — 3)/2. The critical momentum 

is calculated in units of rrieC as Pc = ^'^^iocaili^'^B\ 
where the gyrofrequency is given by vb = eB / {27rmec) 
and B is the magnitude of the magnetic field in the plane 
of the sky. 

A ray casting engine was used to convert volume emis- 
sivities into images with the source set at a user-defined 
distance. In the engine, rays were cast normal to an 
image plane that was rotated to a user-defined orien- 
tation with respect to the grid. Trilinear interpolation 
was used to determine the emissivity at regular inter- 
vals along the ray. The sum of emissivities along a line 
of sight determined the intensity, and the intensity was 
then converted into a fiux by applying the appropriate 
redshift correction and multiplying it by the solid angle of 
a pixel determined by the luminosity distance Dl of the 
source. For both X-ray and radio observations presented 
here, the computational domain was set at a redshift of 
z = 0.0594 {Dl = 240 Mpc) . which is approx imately the 
same as the Hydra Cluster (jWise et a/.ll2007f ). 

3. CLUSTER PROPERTIES: INITIAL CONDITIONS 

To set the stage for analysis we outline now and illus- 
trate in Figures [3] -[5] the initial properties of the ICM for 
our cluster before the jets were launched. Figure [3] shows 
from the simulation initial conditions azimuthally aver- 
aged radial profiles for gas density and pressure (along 
with the derived temperature and sound speed), as well 
as the ICM fiow speed with respect to the DM core and 
the magnetic field strength. The variance and extremes 
of each quantity are also shown. The cluster has a dense 
core with a central ICM density of p/cM,o = 1.05 x 10~^^ 
g cm~^ and pressure P/cM,o = 2.69 x 10~^^ dyne cm~^. 
This gives a core temperature of Ticm,o = 1.6 keV, for 
the assumed mean molecular weight of /i = 0.6, and a 
sound speed of cicm,o = 653 km s~^. The azimuthally 
averaged temperature falls off slowl y with radius, whic h 
is typical for non-cool core clusters (jHudson et a/.|[2QTQ[ ). 
This is not surprising, as radiative cooling was not in- 
cluded in the GADGET simulation that evolved the clus- 
ter to this time. The average magnetic field magnitude in 
the core is Bicm,o = 4.2 /iG, which corresponds to a rel- 
atively high plasma /3/cM,o = 87rP/cM,o/P|cM,o = ^82. 
Figure |4] shows a radial profile of the ICM entropy for 
the initial conditions. This cluster did have a low en- 
tropy core, which is a standard property of SPH sim- 
ulations like the on e that produced the cluster {e.g., 
iMitcheh et"a/1[2009h . 

Although this cluster was selected especially because 
its last major merger was several gigayears in the past 



and at first appearance it seems relatively relaxed, closer 
inspection reveals significant ICM "weather" character- 
istic of sloshing behaviors, thus, reminding us that g676 
has continued to interact with smaller DM halos and to 
accumulate mass slowly. Note first from Figure [3] that the 
mean magnitude of the ICM velocity in the inner core is 
viCM,o = 142 km s~^ (Mach 0.2), and that it increases to 
^ 400 km s~^ by 100 kpc from the cluster center. For a 
typical sound speed of 450 km s~^ outside 100 kpc, these 
motions approach Mach 1. They are primarily relatively 
lar^e scale fiows across the cluster, since by contrast, 
IZhuravleva et all ()2011f ) report isotropic, random veloci- 
ties with respect to the mean of only ~ 50 km s~^ in the 
core of this cluster, increasing to around ~ 100 km s~^ 
outside a 100 kpc radius. The turbulent pressure in the 
core region {Pturb ~ 3 x 10~^^ dyne cm~^), therefore, 
is less than 1% of the core thermal pressure (Figure [3] 
shows profiles of density and thermal pressure). 

Further indications of the sloshing behavior are found 
in Figure [3 which shows orthogonal slices of gas entropy, 
S = P/p^. The initial cluster conditions were rotated 
such that the jet axis from Rl (see Q is in two of the 
planes on the left panel, and the panel on the right was 
created in the same fashion for R2 (see 21 • An inverse 
colormap was used, where lighter colors correspond to 
lower entropy. The jet axis is shown as a dark line from 
the cluster center. Overlaid are ICM velocity vectors 
whose length corresponds to the gas speed. Both views 
show comples fiow patterns and highly asymmetric en- 
trop y distributions, whi ch are classic signatures of slosh- 
ing (jZuhone et a/.ll2010f ). 

Despite these fiow patterns this cluster does still qual- 
ify by common observational measures as an example of 
a relaxed system. A synthetic X-ray observation of the 
initial cluster in the 0.5-8 keV band is shown in Figure [6] 
along a line of sight roughly orthogonal to the mean ICM 
velocity in the cluster core. This line of sight would also 
be along the jet axis from the Rl simulation (see 21 • ^1" 
though there is some asymmetry of the cluster core for 
this particular line of sight, other observational orienta- 
tions lead to more symmetric images. As a quantitative 
measure we construct from this image the concentration 
parameter, c, sometimes used to separate relaxed clusters 
from those that have recently undergone mergers {e.g., 
.Santos et a/. 2008): 



S{R < lOOkpc) 
S{R < bOOkpc) ' 



(9) 



where S is the X-ray fiux in the aperture. A relaxed 
cluster would have a more concentrated core, so a higher 
value of c , than a cluster that was recently disturbed in 
a merger. Cass ano et all (|2010') find that a value c > 0.2 
separates clusters without radio halos (not recently dis- 
turbed) from those with radio halos and a recent major 
dynamical disturbance. From the X-ray observation in 
Figure [6l the initial cluster conditions used for this study 
have a concentration parameter of c = 0.66, so the X-ray 
profile would not be considered dynamically disturbed 
by this observational method. This measurement was 
repeated for two other orthogonal lines of sight through 
the cluster with variation in c of less than 1%. For refer- 
ence we note that it has been pointed out ([Cassano et all 
|2010[ ) that the concentration parameter is less sensi- 
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Figure 3. Radial profile of ICM properties from the initial conditions. The mean with standard deviations are the solid lines surrounded 

by the colored band, and the minimum and maximum values at each radial sample are the dashed lines. 

Note: the lower variation values for the magnetic field magnitude were consistent with zero for some radii and were then truncated. 



tive to projection effects than the so-cahed power ra- 
tio based on multipol e expansions of the X-ray bright- 
ness di stribution fe.a..lBuQte fc Ts ai 1995; Jeltema e^ al\ 
2QQ5I: IVentimJgha et al\ I2QQ8I : lBQhringer"er"a/J I2Q1QI : 
Cassano et a/.|[2QlQf ). 
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Figure 4. Radial profile of ICM entropy from the initial conditions. The mean with standard deviations are the solid line and color band, 
and the minimum and maximum values at each radial sample are the dashed lines. 




Figure 5. Slices of gas entropy, S = P/ p^ , in three orthogonal planes from the initial conditions. The panel on the left shows the initial 
conditions rotated such that the jet axis from Rl, seen as a dark black line from cluster center, is at the intersection of two of the planes. 
The panel on the right is the same for R2. Light to dark colors correspond to low to high values of gas entropy. Overlaid in both panels 
are velocity vectors in the same three planes as entropy. Longer vectors correspond to higher speeds. 
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Figure 6. Synthetic X-ray observation from 0.5 to 8 keV of the 
initial cluster conditions. The image is scaled logarithmically over 
^ 3 orders of magnitude. The line of sight for this obser vatio n was 
along the jet axis from the Rl simulation described in ^4.11 

4. JET SIMULATION RESULTS 

Two AGN outflow simulations for the parameters de- 
fined above in §2.21 were run with the jet cylinder at dif- 
ferent orientations on the grid. Both AGNs were place 
at rest in the cluster center in order to distinguish the 
influences of the dynamic ICM from host galaxy motions 
within the cluster. The first simulation, Rl, had the jet 
cy Under orientation chosen arbitrarily, but so that the jet 
axis was not along a primary grid axis. As it turned out 
the jets were directed roughly orthogonal to the mean 
ICM velocity inside a 200-300 kpc cluster radius, which 
represents the final scale of the cavities formed in the 
ICM by the AGN outflows. The second simulation, R2, 
had the jet cylinder in the orthogonal plane to Rl and 
in a direction more nearly aligned with the mean interior 
ICM flows. The total elapsed time for each simulation 
was 200 Myr. The jets cycled on and off with a 50% 
duty cycle for 6 periods of 26.2 Myr per cycle and then 
left off after 160.5 Myr of this activity, allowing the AGN 
inflated bubbles to evolve without more energy injection 
for about an additional 40 Myr. We note in passing that 
a steady jet of the same properties over this same time in- 
terval would propagate off the grid. A steady jet would 
create a structure roughly twice as long as those pro- 
duced in the two simulations we present. Additionally, 
the X-ray cavity produced by a steady jet would be sig- 
nificantly less round in shape as compared to those pro- 
duced by the intermittent jets in the simulations shown 
here. 

4.1. Simulation 1 (BA) Evolutionary Summary 

Figure [71 shows color, Cj (= jet mass fraction), volume 
renderings from the Rl simulation at five representa- 
tive times (t = 34.4, 49.2, 98.4, 147.6 and 196.8 Myr). 
The views are orthogonal to the jet axis. All the panels 
project the same view except the lower right, as indi- 
cated by the stick diagrams in the two bottom panels. 
A line segment in the lower right panel indicates a pro- 
jected 100 kpc length. To aid structural registration the 
location of the jet cylinder is shown independent of jet 



activity. 

At 34.4 Myr (upper left), the AGN is active for the sec- 
ond time in the simulation, so the jets are collimated, and 
there are small lobes formed from the first active AGN 
phase. As mentioned, these jets drove through a strong 
cross flow in the ICM associated with sloshing motions. 
Already at this time the lobes are obviously asymmetric 
and no longer aligned with the jet axis. At 49.2 Myr 
the AGN is inactive, and it is apparent that the AGN 
plasma is concentrated at the ends of the lobes. The 
panels at 98.4 Myr and 147.6 Myr both show times when 
the jets are on. In both cases the visible jets clearly 
are not as well collimated as they were at 34.4 Myr. The 
ICM pressure near the launching region is approximately 
a factor of two lower at these times than the value at 34.4 
Myr. Recall from §2.21 that the jets were launched with 
a fixed pressure related to the initial ICM, so the jets 
are now over-pressured as they emerge, causing them 
to expand. Since the internal jet Mach number is low 
(Mj ^ 1.2), re-collimation is slow. In addition, bow 
shocks have formed ahead of the new outflows, as well as 
jet termination shocks, further enhancing the dispersion 
of jet plasma into the cavity associated with the previ- 
ous activity. The asymmetry between the lobes and their 
misalignment with the jet outflows are quite obvious, es- 
pecially in the final views at 196.8 My, sr. The lobes 
extend roughly the same distance from the cluster center 
(300 kpc), but there are different offsets of the lobes from 
the jet axis. The most significant differences between the 
two lobes are the deflection in roughly the —z" direction, 
so that viewed roughly along the y" axis they present a 
"C-shaped" morphology, and from another, orthogonal 
perspective, a flattened shape of the left lobe compared 
to the right lobe, as seen in the bottom right panel of 
Figure [71 Finally, we note that the AGN plasma is very 
nonuniformally distributed inside the lobes, with some 
significant mixing of ICM and AGN plasma in several 
regions, made evident by the gray (turquoise for color fig- 
ures) regions in all the figure panels. Such mixing should 
have various observational consequences, including the 
Faraday rotation behaviors along sight lines through the 
lobes, as we will discuss in a subsequent publication. 

To that point, however, we mention here that magnetic 
fields, while not dynamically dominant in these simula- 
tions, are essential in the production and transfer of radio 
emissions. In those roles they can help to illuminate dy- 
namical properties of the AGN/ICM interactions. Figure 
m shows volume renderings from the Rl simulation of the 
color, Cj, along with selected magnetic field lines colored 
dark to light by magnitude during an active AGN phase 
at 137.8 Myr. In the left panel the jet to the right is com- 
ing towards the observer and the left jet is directed away. 
From this perspective the AGN outflows are deflected 
away from the observer. Also in the left panel, ICM field 
lines can be seen stretched and dragged outward with 
AGN plasma into the leftward lobe formed from previous 
activity. Those ICM field are actually draped around the 
AGN plasma. The field lines colored white show where 
the stretching from the jet outflow has most intensified 
the field. OJIO reported a similar result for their inter- 
mittent jet simulations. There is a stray magnetic field 
line seen in the left panel of Figure [51 Despite attempting 
to isolate field within the jet, this line is rendered due to 
the blending of the jet and ICM magnetic fields in the 
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Figure 7. Volume rendering of the color variable, Cj, from Rl. For color the color version of this figure, the color map has ICM plasma 
in red, equally mixed ICM and AGN plasma in turquoise, and predominantly AGN plasma in yellow to white. The grayscale version of 
this figure shows AGN plasma in white and ICM plasma in black. 
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collar region around the jet cylinder. In the right panel 
the magnetic field inside this same lobe is shown from 
a slightly different perspective wrapping around regions 
consisting primarily of AGN plasma. Field lines are also 
evident tracing the outline of the lobe. 

4.2. Simulation 2 (BJl) Evolutionary Summary 

The evolution of the R2 simulation is shown in Fig- 
ure [HI again through volume renderings of the color Cj . 
Many of the same morphological characteristics of Rl 
visible in Figure [7] can be seen in this rendering of R2, 
as well. Just as in the Rl case the R2 jets appear well 
collimated early in the simulation, but again they over- 
expand as the propagate through the lower pressure cav- 
ities. The AGN plasma in the lobes has a qualitatively 
similar distribution to Rl. Again the AGN plasma is 
primarily deposited near the outer edges of the lobes. 
The most striking difference between the Rl and R2 
morphologies is that while the lobes in the Rl simula- 
tion extended about the same distance from the AGN, 
the lengths of the two lobes in R2 are of very differ- 
ent length. The lobe on the right side of Figure [9] does 
not extend nearly as far into the cluster as the lobe on 
the left side in these images. This is evident even in 
the early, 34.4 Myr image. The asymmetry continues to 
increase throughout the simulation. At 196.8 Myr the 
difference in extent is approximately 100 kpc or roughly 
50%. Recall that these jets are directed in an orthogonal 
direction to those in the Rl simulation. In fact, the jet 
axis is roughly aligned parallel and anti-parallel to the 
prevailing ICM velocity in the inner cluster mentioned 
previously. The relationship between cavity morphology 
and cluster weather will be discussed in more detail in 

The structure of several magnetic field lines during jet 
activity at 137.8 Myr from R2 are seen in Figure [TOl 
Similar to the structures in Figure [HI the toroidal jet field 
is shown in the left panel, along with ICM field lines that 
have been stretched and wrapped over AGN plasma from 
the last period of activity. Magnetic field lines inside the 
leftward lobe are visible in the right panel. Those lines 
also follow the outline of lobe structures, as they did in 
the Rl simulation. The fields seen here appear to be less 
tangled than those seen in Figure [H This same lobe is 
the one on the left-hand side of the images in Figure [9) 
Its distortions are not as extreme as the deflected jet and 
lobe shown from Rl in Figure [H 

4.3. Synthetic X-ray Observations 

Figure [TT] shows 0.5-8 keV X-ray observations of Rl 
at various epochs taken along a line of sight that was 
perpendicular to the jet axis. The resolution is 1 arc sec 
per pixel at the assumed 240 Mpc distance, and each 
obse rvation was divided by the best-fit double /3-profile 
(seelMendvgral et al\\2^1\\ \ to emphasize the X-ray cavi- 
ties. In these images, the lower X-ray cavity is associated 
with the left-hand lobe in Figure [3 The observation of 
the initial conditions shows significant departures in sur- 
face brightness from the double /3-profile, particularity 
at large distances from the cluster center (see, also Fig- 
ures [5] and [6|) reminding us, once again that the ICM is 
not truly relaxed and the cluster is not isolated. At 65.6 
Myr, after two full AGN periods, the jets are active and 



and the inner pair and outer pair of cavities are seen sep- 
arated by a thin bright rim. These featu res are similar 
to th e multi-cavity system in Abell 2052 (iBlan ton et al\ 
120091 ). The observation at 131.2 Myr shows the same 
morphology with the pairs of cavities on larger scales. 
Five complete AGN cycles have been executed by this 
time. Ripples are seen emanating from the outer pair 
of cavities, which are outlined by bright rims that are 
~ 30% hotter than the surrounding gas. The ripples 
were seen for th e intermittent simulation presented in 
Mendvgral et al\ ([2011), and they are similar to t he rip - 
ples observed in the Perseus Cluster (jFabian et a/.ll2003l ). 
Additionally, bright rims hotter than their surroundings 
have been observed i n Centaurus A (iKraft et al\ l2003l 
■200?) and NGC 3801 (|Croston et a/1l2007D . but for other 
objects, th e bright rims are, i n fact, cooler than their sur- 
roundings rtPiehl et a/.l l2008\ 

The outer Rl cavities at 131.2 Myr are distorted and 
no longer aligned with the axis of the inner pair of 
cavities. These distortions are largely due to cluster 
"weather" , as mentioned above and discussed more fully 
in ^ The approximately Mach 2 bow shock is well sep- 
arated from the cavities at this time. By the time of the 
final observation, at 196.8 Myr, the AGN has been inac- 
tive for 40 Myr, and there is no longer an apparent pair of 
inner cavities. The cavity system appears to include only 
one pair that extend ^ 300 kpc from the cluster center 
and are ~ 150 kpc wide. The cavities are still outlined 
by bright rims, and the ripples can still be seen. The 
lower cavity shows the most distinct distortions. The 
bow shock, now Mach ~1.3, is still visible. 

Figure [12] shows 0.5-8 keV observations of the R2 simu- 
lation, also divided by best-fit double /3-profiles. In these 
images the fine of sight is again perpendicular to the jet 
axis. Many of the same feature seen in the Rl observa- 
tions are present, such as the ripples from the intermit- 
tent AGN activity, bright rims around the large cavities 
and a bow shock with a similar Mach number to the Rl 
case. A unique and very interesting feature of the R2 ob- 
servations, however, is that the 65.6 Myr and 131.2 Myr 
observations show an incomplete pair of inner cavities 
whereas Rl showed a complete pair. The rim around 
the upper-inner cavity is visible in these two images, but 
the lower-inner cavity is minimally outlined with a faint 
rim. It is unlikely that this cavity would be seen in a 
real observation of such an object with finite sensitivity, 
and only three of the four cavities might be identified. 
In the last image at 196.8 Myr the inner cavities are ab- 
sent, and only a single pair of cavities with bright rims, 
that are less noticeable than those from Rl, are visible. 
The lower cavity is significantly shorter and more circu- 
lar in appearance than the upper cavity. As mentioned 
previously and discussed in more detail in ^ this asym- 
metry also results from large scale ICM fiows (weather) 
present in the cluster initial conditions. As mentioned 
previously the R2 jets are more or less aligned with this 
fiow, whereas the Rl jets were roughly orthogonal to the 
mean interior ICM fiow. 

4.4. Synthetic Radio Observations 

Figures[T3land[T4lshow at 131.3 Myr and 198.2 Myr the 
178 MHz synthetic synchrotron intensity observations of 
the Rl and R2 simulations, respectively. These obser- 
vations are of the same line of sight as the respective 
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Figure 8. Volume rendering of passive tracer Cj and selected magnetic field lines from Rl during the active AGN phase at 137.8 Myr. 
Field lines are colored by magnitude with brighter field line colors representing stronger fields, left: A view showing the tightly coiled jet 
field near the launching region. ICM field lines are seen stretched over the AGN plasma from a previous AGN cycle, right: Field lines in 
this view within the leftward jet lobe in the other panel are seen following the outline of the lobe. 



X-ray observations discussed above, but they are of from 
a smaller volume to show more detail. We note to avoid 
confusion that in all the images the jet launch cylinder 
remains artificially illuminated. The left-hand panels in 
both figures were configured with 1.7 arc sec beams (pix- 
els) at the assumed 240 Mpc distance, and the right-hand 
panels are the same observations but Gaussian smoothed 
over 6 pixels to simulate a 10 arc sec beam. The line of 
sight through the grid in Figure [13] (Figure [14]) is the 
same as the X-ray image in Figure [TT] (Figure [T2]) . Such 
low frequency radio observations are of particular inter- 
est, because they highlight low energy CRs that provide 
a reasonable tracer for AGN plasma, since their losses 
are not large. For a typical field value in the lobes of ^ 
1.5 /iG much of the emission at 178 MHz is from ~ 3 
GeV electrons, whose radiative lifetimes against inverse 
Compton and synchrotron losses in this context are ~ 
200 Myr, assuming a small cluster redshift. 

For both simulations at 131.3 Myr, the bright ac- 
tive jets are seen in the figures penetrating their asso- 
ciated lobes. Emission from previous AGN activity ex- 
tends further out, but is considerably dimmed because 
of radiative aging of the CR population. These images 
closely resemble a ^^double- double" radio galaxy (DDRG) 
([Schoenmakers et a/.ll200Q) such as 3C 293, which also 
show s a misalignment between the inner and outer dou- 
bles (|Joshi et al\\2011\ ). At the end of the simulations, 
most of the DDRG resemblance is gone because roughly 
50 Myr has passed since the AGN was last active. Rl 
shows a distinctly bent morphology, similar to a WAT 
radio source. The R2 observation in Figure [14] has a sec- 
ond pocket of emission in the lower lobe, which has not 
extended as far from the cluster center as the upper lobe. 
The radio emission shown in Figures [T3] and [T4] was con- 
tained within the boundaries of the X-ray cavities seen 
in Figures [11] and [11] respectively. 

In addition to qualitative morphological evaluations 
that reveal simple dynamical relationships, the above 



synthetic observations can be examined quantitatively 
to go od effect, as well. For example, iTregillis et al\ 
(|2002l ) used such observations at multiple frequencies to 
test observational methods for measuring magnetic field 
strengths and properties of the CR population. In ad- 
dition to total intensity maps, our synthetic observation 
tool can also take advantage of the vector magnetic field 
and CR electron plasma information in the simulations 
to construct polarization maps that include Faraday ro- 
tation. These can be used, for example, to test tech- 
niques for measuring the structure of magnetic fields in 
the AGN and in the ICM surrounding it. We defer those 
studies to subsequent publications. 

5. THE ROLE OF CLUSTER "WEATHER" 

The asymmetries between the jet lobe and X-ray cav- 
ity pairs in these simulations clearly demonstrate that 
asymmetries in the properties of the ICM, that is, clus- 
ter "weather patterns", can strongly influence the mor- 
phology of radio lobes and ICM cavities. Since most ra- 
dio galaxies in cluster environments show some degree of 
distortion, it is useful to understand how to establish the 
dynamical links to observed AGN outflow distortions as 
metrics of the ICM structure. In our jet simulations we 
can exclude direct gravitational forces as contributors to 
nonradial distortions of the AGN plasma and entrained 
ICM, since the gravitational potential in action during 
the simulations was spherically symmetric. Nonspheri- 
cal gravity was essential, of course, in the development 
of the ICM conditions present at the start of our simu- 
lations, especially those we describe as sloshing. 

A simple model can illuminate the main character of 
the AGN/ICM-weather interactions at work during our 
simulations. For this purpose we can picture the AGN 
outflows as being ejected into a wind at some angle to 
the jet axis. It is suflicient for our purposes to look at 
the two limiting cases having the winds aligned with the 
jet axis (i.e., "head" and "tail" winds) and orthogonal 
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49.2 Myr 




Figure 9. Same as Figure [71 for R2. 

to the jets ("cross" winds). For simplicity we assume a 
uniform ICM and steady jets. In doing this we should 
remember that our simulated environment is nonuniform 
and the simulated jets are intermittent. Still, the model 
can provide some insights and consistency checks. 

In the head and tail (aligned) wind case that applies 
approximately to our R2 simulation we expect the length 
of the jet structure, ^, on the tail wind side {v^ > as- 
suming Vj > 0) to be extended in comparison to a sta- 



tionary ICM, while the head wind side {v^ < 0) will 
be shortened. A modified version of the classic calcula- 
tion to balance the jet thrust agains t ram pressure from 
the ICM as the jet advances (Blan dford fc ReesI 119741 ) 
provides a simple estimate of the wind's infiuence. The 
model estimates, if the local ICM density is picM and 
the jet thrust is Fj, that the rate at which the jet "head" 
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Figure 10. Similar to Figure [8] for R2. 



extends, Vh > 0^ is 



Vh ^ VhO 



= ^^0 1 



K\ V Fj 



-Vx 



1/2- 



(10) 



where r/j is the cross sectional radius of the head, and 



VhO 



i + Vx' 



(11) 



with X = {"^h/TjY^picM/ Pj)' In the apphcable,"hght 
jet" hmit, y/x >> 1 this gives the intuitively obvious 
result, 



Vh ^ VhQ + v^. 



(12) 



Setting the (observed) ratio of lengths, 5 = i-\-/i- ~ 
{vho + I'^wD/ivho -\vw\), one can find, noting that Vho is 
the average extension rate on the two sides, an estimate 
of the wind speed as, \vw\ ^ Vho{^ + l)/(^ ~ !)• 

At the end of our R2 simulation {t ~ 200 Myr) one 
jet cavity extends roughly 300 kpc from the AGN, while 
the other reaches a distance only about 200 kpc from the 
AGN (see Figure [9|). Thus, we can estimate S ~ 1.5. 
The average length is then about 250 kpc, which over 
the approximately 200 Myr life of the AGN activity, cor- 
responds to an average extension rate (~ Vho) around 
2500 km s~^ (approximately Mach 5). Using the rela- 
tion derived in the previous paragraph this would pro- 
duce an estimate for the ICM wind speed, Vuj ~ 500 km 
s~^. On these scales the actual ICM fiow speed is around 
400 km s~ (see Figure [3]), so given the simplicity of the 
estimate, consistent with this result. 

The orthogonal, cross wind case, related to our Rl 
simulation, is familiar in NAT and WAT contexts. In 
particular a transverse ram pressu re from the wind will 
defiect the jet and its cocoon (e.g.. Begelman et a/.lll979l : 
[jones fc Owenl 119791) . In the absence of a wind the jets 
will deposit momentum, energy (including pressure) and 



AGN plasma in the head region. The energy and plasma 
will fiow backwards and expand in response to the added 
pressure to form a AGN plasma cocoon, which is rep- 
resented by X-ray cavities in clusters. The cavities are 
reasonably clo se to pressure balance with the surround- 
ing ICM fe.g. JMcNamara fc Nuls"enl 120071 OJIO), so the 
rate of lateral expansion will be of the order of the ICM 
sound speed, cjcm- The Rl cocoons at the end of the 
simulation have, on average, transverse radii. Re ~ 100 
kpc, corresponding to mean expansion rates around 500 
km s~^. That is close to cjcm and consistent with this 
expectation. 

In a cross wind the cocoon plasma will experience a 
ram-pressure-based lateral force, F^ ~ 2picMv'^^Rc^ 
where i is the length of the cocoon and Re is its trans- 
verse radius. The mass in the cocoon is roughly peiirRl, 
with pe the density in the cocoon. So, the time for F^ to 
accelerate the cocoon to the wind speed is roughly, t^c ^ 
{7r/2){pe/picM){Me/M^) te, whcrc M^ = Vyi/cicM and 
te ^ Rc/{MeCicM) IS the time for the cocoon to ex- 
pand. With both Me ~ 1 and M^ ~ 1 the cocoon 
will then reach terminal speed at rest in the ICM on 
a timescale less than it takes the coccon to expand. On 
the other hand, if the wind speed exceeds the expansion 
speed; namely, if M^) > Me, the cocoon will be stripped 
from the jet, leaving a "naked" jet. Our Rl jets do 
possess cocoons, since M^ < Me . Assuming that the 
cocoons move across the jet path at the wind speed, we 
can estimate the down wind displacement angle of the 
cocoon compared to the jet direction at its source, 6d, as 
tan^(i ~ Vw/vh- For the Rl structures at the end of the 
simulation, this gives tan6>(i ~ 0.15 — 0.2, so ^^i '^ 10 — 15 
degrees. That is consistent with the deformations visible 
on the left side of Figure [71 

For completeness we recall that a naked jet is defiected 
by a cross wind ove r a characteristic l ength, £d^ (e.g. 
iBegelman et a/.lll979l : rPorter et a/.ll2009f ). 



id{naked) 



Pj 



Mf 



PICM V. 



2 'J 



ICM 



Ml 



(13) 



16 




0.06 



0.24 



0.54 



0.96 



1.5 



2.2 



2.9 



3.8 



4.9 



Figure 11. Synthetic X-ray observations from 0.5 to 8 keV from four epochs of Rl. Each observation was divided by a best-fit /3-profile 
to emphasize the X-ray cavities. 



The expression on the far right emphasizes that the de- 
flection length, id of a naked jet in pressure equilibrium 
with its environment depends only on the jet radius and 
the relative Mach number in the jet flow to the Mach 
number of the wind. If the jet is cloaked it will be de- 
flected by a transverse pressure difference across the co- 
coon resulting from F^. T he deflection length in this 
case can be estimated to be (j Jones fc Owen|[l979f ). 



£d{cloaked) 



Pj 



PICM V^ 



-Rr 



Re 



id{naked), (14) 



so the bending will be more gradual for a given wind. 
In our case the expected ratio is roughly an order of 
magnitude, consistent with the simulation results. 



6. CONCLUSIONS 

We have presented two MHD simulations of AGN jets 
inside of a galaxy cluster that was adapted from an SPH 
cosmological simulation to explore the effects of cluster 
"weather" on the jets and lobes. The important results 
from this work are: 

• Evidence that ICM flows in these simulations were 
sufficient to produce asymmetries between the op- 
posing AGN jets and lobes was presented. Al- 
though pressure variations in the ICM (unrelated 
to the presence of a gravitational potential) and 
magnetic stresses were present, they were either 
not as large as the variations in the ICM motions 
or were very localized. The existence of significant 
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Figure 12. Same as Figure [TT] for R2. 



0.96 



"weather" in the cluster used for these simulations 
is interesting because it was specifically selected as 
a relatively relaxed system. ^ showed that typical 
flow speeds in the cluster were one the order of 400 
km s~^, which appear to be sufficent to modify the 
jets and lobes. This result is in agreem ent with 
the predictions of lHardcastle et al\ (J2005l ). who ar- 
gued that flows as low as 100-300 km s~^ could cre- 
ate bent WAT-like morphologies for jets similar to 
those in our model (high velocity and low density). 
These results also show that the AGN outflows act 
as probes of bulk motions in the ICM, even in a 
seemingly relaxed cluster. Simple dynamical mod- 
els utilizing the influence of the ICM ram pressure 
on the AGN cocoons seem to give useful estimates 



of the ICM properties. This makes it possible to in- 
fer velocity structure information in environments 
where direct measurement is not possible. 



• Magnetic field lines anchored in the ICM are 
stretched outward by the intermittent jets. This 
causes field lines to wrap around the jets and AGN 
plasma in the lobes, and in localized regions, the 
stresses associated with the magnetic field were on 
the order or even greater than the inertial stresses. 
Such regions serve as a reminder that even for the 
"high-/3" plasmas modeled here. Maxwell stresses 
can significantly effect the evolution of a plasma 
and emphasize the importance of modeling these 
systems with MHD. 



18 



131.3 Myr 

T 


1 


T 


J 


196.8 Myr 

i 

1 


1 

• 

1 00 kpc 


^^^^^^^^^^^^^^H 



4,66e-07 



3,27e-06 



l,45e-05 



5,87e-05 



2.36e-04 



Figure 13. Synthetic observations of synchrotron emission from Rl at 178 MHz. The units on the colorbar are Jy beam 



dth a 1.7 



arcsec beam (pixel), left: Full resolution observation, right: Same observation smoothed with a Gaussian kernel over 6 beams (pixels). 



• Synthetic X-ray observations of these simulations 
show complex cavity structures associated with the 
intermittent jets and ICM motions. Multiple cavity 
pairs were "observed" as restarted jets pushed into 
cavities formed by previous outbursts. In one ex- 
ample given, a new cavity forming inside a larger 
cavity did not have an obvious companion cavity 
from the opposing jet. Bright rims were seen out- 
lining the cavities, and they were associated with 
the injection of energy by restarting jets. Addi- 
tionally, cavities were seen misaligned with the jet 
launching region due to the deflection of the jets 
by ICM flows. 

• Synthetic radio observations made during periods 



of jet activity resembled "double-double" radio 
galaxy sources with a bright pair of jets interior 
to fainter sources further out in the cluster. This 
state was temporary, however, as the cosmic ray 
electrons from the new outburst would eventually 
reach the end of the outer lobe. The ICM motions 
created a WAT morphology of the lobes in one sim- 
ulation. 
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